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Abstract 

We describe a method of solving quantum field theories using operator 
techniques based on the expansion of interacting fields in terms of asymptotic 
fields. For bound states, we introduce an asymptotic field for each (stable) 
bound state. We choose the nonrelativistic hydrogen atom as an example to 
illustrate the method. Future work will apply this N-quantum approach to 
relativistic theories that include bound states in motion. 



1 Introduction 



We describe a method for solving of quantum field theories (QFT) that has wide 
applicability, and is analogous to methods developed to solve problems in classical 
theories. To find approximate solutions of a theories with a given Lagrangian, the 
method approximates the Lagrangian fields by a series of normal-ordered asymptotic 
fields that terminate with a product of up to asymptotic fields. This is the 
"N" of the N-quantum approximation. R. Haag introduced this expansion pQ. We 
modified his work by introducing asymptotic fields for the stable bound states, if 
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any, in the theory. We derive analogs of the Schrodinger equation that have solutions 
when the theory has bound states. 

We apply this method to the hydrogen atom to illustrate how this works for 
bound states. Following C.N. Yang and D. Feldman [2], we use the QFT analog 
of the solution of classical equations of motion in which the the classical equations 
govern a set of functions while the QFT equations of motion govern a set of operators. 
As stated above we take account of the operator aspect of quantum field theory by 
using a complete and irreducible set of asymptotic fields, either the in or the out 
fields [Tj. A complete set of fields is one that creates all the states of the Hilbert space 
by acting repeatedly on a given vector, often the vacuum, called the cyclic vector. 
(More accurately, the complete set needs only to approximate any vector when 
polynomials in the operators, integrated ("smeared") using smooth functions, act 
on the cyclic vector.) An irreducible set of fields is a set that has no proper invariant 
subspace. By Schur's lemma, any operator that commutes with all operators in an 
irreducible set must be a multiple of the identity operator. Saying that the in (or 
out) fields comprise a complete (and irreducible) set, means that any operator in 
the theory can be expressed in terms of them. 

The in (or out) fields have free-field commutators, obey free equations of mo- 
tion, and the different in (or out) fields commute or anticommute with each other 
everywhere in spacetime. Each of these sets of asymptotic fields is completely known 
once the masses, spins and quantum numbers of the fields in a given set are defined. 
Thus, either set serves as a collection of standard building blocks from which we can 
construct solutions of the operator equations of motion. 

To solve the equations of motion using the Yang- Feldman approach, we expand 
the fields that appear in the Hamiltonian or Lagrangian in normal-ordered series 
of in (or out) fields. (This is the Haag expansion.) To determine the c- number 
amplitudes (the Haag amplitudes) that are the coefficients of the normal-ordered 
terms, we insert the expansion in the operator equations of motion, again normal- 
order the asymptotic fields, and equate the coefficients of corresponding (linearly 
independent) normal-ordered terms. 

In this paper, we illustrate this technique for the case of the nonrelativistic 
hydrogen atom [3j. The amplitudes for the bound states satisfy the Schrodinger 
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equation for the hydrogen atom and, of course, we find the usual spectrum of bound 
states. This shows that the relevant Haag amplitudes are the wave functions of the 
bound states. 

Our method, which is based on the Haag expansion, is entirely independent 
of the Bethe-Salpeter equation. In contrast to the Bethe-Salpeter approach, our 
method contains no spurious solutions and no negative norm amplitudes. The 
method can be used for bound states in relativistic theories, despite that the ampli- 
tudes depend on the same number of kinematic variables as appear in nonrelativistic 
wave functions. There are no relative time intervals in the relativistic version of this 
method. 

We applied this formulation previously to a model of the deuteron [1] in an 
approximation where this method is similar to that of the Gross equation [5], but 
differs from it in higher-order approximations. 

In this paper, we neglect two related issues that are relevant for the Coulomb 
interaction. One is that, because the Coulomb interaction is long ranged, spatially 
separated particles never move independently [0]. This does affect any questions 
addressed in our study. The other point is that the Coulomb field of charged particles 
is not attached to the particles, as suggested by P. A.M. Dirac [7]. Nevertheless, this 
too does not appear to cause any difficulty for the conclusions drawn from this work. 

Although our method does not require expressions for interpolating (Heisen- 
berg) in and out asymptotic fields for bound states, we show how to construct such 
interpolating fields for the bound states of the theory (in our case, the states of the 
hydrogen atom) in terms of the fundamental fields that enter the Hamiltonian of 
the theory. This construction, which uses a convolution of the bound-state wave 
functions with the constituent fields, contrasts with the constructions of K. Nishi- 
jima |H] and of W. Zimmermann [S], which use products of the fundamental fields 
at the same point. 

We plan to apply the construction of interpolating fields for bound states of 
baryons and mesons to calculations based on quark models within quantum chro- 
mo dynamics. 

We note that when we use in (out) fields as the asymptotic fields, we have 
retarded (advanced) boundary conditions on the amplitudes (Green's functions) that 
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provide dispersion relations for these amplitudes, which G. Kallen [10] discussed 
in the context of perturbation theory for Yang-Feldman equations with retarded 
boundary conditions. 



2 Asymptotic fields and the Haag expansion 

Asymptotic fields have been part of quantum field theory at least since the work of 
H. Lehmann, K. Symanzik and W. Zimmermann [11]; nevertheless, there are still 
misconceptions that should be cleared up. Asymptotic fields are free fields in theories 
that have no long-range potentials, massless particles or confinement, because, in 
such theories, we describe particles either by: (i) separated asymptotic fields (at large 
time scales) or (ii) as a bound state, also represented by an independent asymptotic 
field. In either case, the exact eigenstates can be labeled by the quantum numbers 
of the free particles. The (weak) limits for t — ?■ ±oo are the out or in fields that 
construct the eigenstates at the corresponding limiting times [12]. The nontrivial 
unitary relation between these fields is given by the ^-operator. 



The asymptotic fields at finite times are the limiting fields propagated to finite times 
according to the free equations of motion. The free-field property of the asymptotic 
fields does not depend on any unphysical "adiabatic switching off" of interactions. 
Even when the potentials are long-range (and there are massless particles in the 
context of the field theory), we can use often use the usual asymptotic fields, as is 
the case for the example considered in this paper. 

The limits that define the asymptotic fields are subtle. The relations that 
appear in some books. 



are ill-defined. We discuss the proper (weak-limit) definition of the asymptotic fields 
in the Appendix. 




(1) 



y4(x, t) A(°"*)^" (x, t), t^±oo 



(2) 
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3 Nonrelativistic model of the hydrogen atom 



The fundamental fields are the spin 1/2 electron, ea{x,t), and the spin 1/2 proton, 
Pa(x, t), fields. The a indices label the spin components of the fields. We assume 
that these fields obey the usual equal-time canonical anticommutation relations, 

[ea(x, t), ej(y, t)]+ = SabS{x - y), (3) 

\Pa{^,t),pl{y,t)]+^5abSiyi-y), (4) 

and that the other equal-time anticommutators vanish. The nonrelativistic Hamil- 
tonian is 

^ = ^ / c?'a;Vxel(x, t) ■ Vxe„(x, t) + ^ J d^xV^pH^, t) ■ VxPa(x, t) 

d^xd^ypliy, t)ei(x, t) _ e„(x, t)pb{y, t), (5) 

l-'^ y I 

where m and M are the electron and proton masses. (We have not included spin- 
dependent terms in the interaction.) The equation of motion for the electron field 
is 

idtea{^,t) = -— V^e„(x,i) - j d^ypl{y,t)pb(y,t)j^^-y^ea{x,t). (6) 
The proton field obeys a similar equation. 

4 Derivation of the Schrodinger equation for the 
hydrogen atom 

We illustrate the N-quantum approach for the hydrogen atom in the approximation 
in which the excited states of hydrogen are stable. We introduce the terms for this 
case into the Haag expansions of the electron and proton fields. The relevant in 
fields (alternatively, the out fields can be used as well) are those for the electron, 
e™(x, t), the proton, pjj"(x, t) and the in fields for each state of the hydrogen atom, 
h^^iy^i. For simplicity, we choose the F = state of the hydrogen atom, so the 
in field, h^j^i,^,, is a scalar boson. (We use the known spectrum of the hydrogen 
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atom. In general, we would use the Schrodinger equation, which we derive below, 
to find any such bound states.) The in fields obey free equations of motion and free 
anticommutation (commutation) relations: For the electron, we have. 



za,er(x,t) = -— Vier(x,t), (7) 

[er (x, t), er^(y, t')]+ = -y,t- t'; 0, m), (8) 

P(x, t; E, m) = — i-^ / dud^k6{uj -E- i^)e-^-*+^k-_ (g) 

The proton in field obeys the analogous equations, with the proton mass replacing 
the electron mass. The in fields for the states of the hydrogen atom obey, 

idtKlmi^.t) = (Enlm ' TTV^'^Xlmi^, t) , (10) 
[^n1m(x, t), /^^"//^/(y, t')]- = Snlm,n'l'm''D{x -y,t- t'] Enlm, Mnlm), (H) 

where Mnim = m + M + Enim is the mass of the hydrogen atom in the state nlm. 
The anticommutators or commutators between different in fields vanish. 

Because the electron is stable the appropriate terms in the Haag expansion 
for the electron field are the electron in field, and for the states of the hydrogen 
atom, the hydrogen in fields. In the annihilation part for the electron, the oper- 
ators that enter the bound-state term are the hydrogen annihilation operator and 
the creation operator for the proton. The quantum numbers balance, in that, for 
example, both e*" and p^^'^K^im annihilate charge — e and lepton number +1. (For 
simplicity, we choose the F = state of the hydrogen atom, so that K'nim represents 
a scalar boson.) As will be shown explicitly below, this term in the Haag amplitudes 
corresponds to the usual bound-state wave functions for the hydrogen atom. We 
choose a scalar hydrogen atom to simplify our discussion. Using the constraint [13] 
from the nonrelativistic boost, k — t- k + v, we find 

e,(x,t) = er(x,t) 

+ E / rmrni^ - y)ea, : pTHj, ^) ^n"n^ ( ! . t) ■■ d\ (12) 

1 J lit ~\ -i- 'J- 

mm 
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The analogous expansion holds for the proton field: 



Pa{y,t) = pT{y,t) 



+ E / - yK, : er(y,t)/.-^( ^^ + j|^y ,t) : d\ (13) 

nlm 

The wave function '^^im matrix element of the Hamiltonian field ea(x, be- 

tween the proton in state and the in states for the hydrogen atom: 

m(P(,(y,t)|ea(x,t)|/i(X,t)),„ = e„b<,^(x - y)5( — / - X). (14) 

m-\- M 

To find the equation for the wave function ipnimi substitute these expansions 
into the equation of motion for the electron field, again normal-order the in fields, 
and examine the terms that contain p'^'^h^nim'- 



•' nlm '''' 
nlm'' l-''^ y\ m-\-IVl 

To find the Schrodinger equation for the hydrogen atom with momentum P, we 
introduce the Fourier transform: 

^n"n^(x, t) = J h^^{P)exp{-z{E^im + ^^)t + iP • ^)d'P. (16) 

In terms of the Fourier transform, the equation is 

Y: [ rmmi^ - y)U^dt + E^i^ + -^) : ^^^(y, t)hZAP) ■ d'y = 

nlm-^ ^^^nlm 



nlm ' 



We evaluate the time derivative on the left hand side using the equation for 
the p^^\ and integrate by parts to get derivatives acting on the wave function. The 
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left hand side becomes 

nlm nlm 

We do not assume the hydrogen atom is at rest; the last term above gives its 
kinetic energy. We are using the kinematics for the nonrelativistic limit of a Lorentz 
covariant theory. (For a Galilean covariant theory, the mass of the hydrogen atom is 
m + M for all states of the hydrogen atom because of the Bargmann superselection 
rule. [2]) The final step is to equate the coefficients of the Pf,"^^^"^ terms. We justify 
this either by noting that these terms are linearly independent of all other types of 
in field products or, more formally, by anticommuting with p^" and commuting with 
h^nim- The result is 



T^^l - 7777^1 - 1 Oi'nlmi^ - y) = (E^lm + 77Tl—)i'nlmi^ " Y)) , 



nlm 



mMe^ 

-t^nlm = ~7T7 [ iTn 9' \^^) 

where we used the known energy levels of the hydrogen atom. We recognize this as 
the Schrodinger equation for the hydrogen atom and identify the il^^im its wave 
functions. 

In terms of the relative coordinate, r = x — y, the Schrodinger equation for 
the hydrogen atom with momentum P is 

' ^ -V^ - i^Wnlmi^) = (Enlm + ^P^)<,„(r), (20) 



fJ'nlm 1 1" I l^nlm 

where l//i = l/m + 1/M and 1/ finim = ^/Mnim — l/(m + M). 

Because 1/finim > 0, the relative wave function for the bound state in motion 
obeys an equation with a (slightly) higher energy than that for the bound state at 
rest. In a Galilean invariant theory, Mnim = M + m for all states of the hydrogen 
atom, and the atom in motion obeys the same equation as the atom at rest because 
of the Bargmann mass superselection rule [1] 
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The bound states of the hydrogen atom, taking into account the spin of the 
electron and the proton, are labeled by nLSJJz, and the N quantum amplitudes 
also have these labels, as well as spin indices a, b for the electron and proton spins. 
For the S = states, the amplitude is 

ra/jSnL S = 0)(r) = e„,V^}/(nL S = 0)(r). (21) 

For the 5* = 1 states the amplitude is 

ra,\inL S = l)(r) = a!s^P}finL S = l)(r), (22) 

where J = L — 1, L, L + 1. The corresponding annihilation operators for the 5 = 
hydrogen atom are 

hi (nLS, J = L;R) = h^^'jU. i^) (23) 

and for = 1, 

hiiriLS, J = L - 1, L, L + 1;K) = 

E v/2JTT(-)^-^-^^ f f / ^ M (24) 

\ J^z Jz ~ i-'z Jz J 

5 Normalization of the wave functions 

The asymptotic fields diagonalize conserved observables such as the Hamiltonian, 
the momentum operators and various charges. Thus, any such conserved quantity 
has the form of a sum of free operators, O/ree{0m}, for the contribution of each 
in field, including those for bound states, to the observable O. For the ip^im^ with 
the electron off-shell, we choose the lepton number L, whose term in the in field 
expansion provides the normalization condition, to be the observable. The lepton 
number in this model is 

L = j t/W(x, 0)e(x, 0) = Lo{e-} + ^ L,{KiJ + ■■■. (25) 

nlm 

We substitute the Haag expansion for e(x, t) given in Eq. flT^ into L, look for the 
terms bilinear in and equate these to Lq{K'^i^} = J d? Rh^^i^K^i^^. We find 

/ E <L(x - y)Cl.(R, 0)p™(y, 0)d^ 

nlm 
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X E ^:nU^-y')rHy',o)h:i,m'i^',o)d'y' 



n'l'm' 



nlm 



(26) 



with R = {my. + My)/(m + M)and R' = {my. + My')/(m + M). The operators 
on the left-hand side contract, (p*"(y, 0)p*""'^(y', 0)) = 5{y — y'). When we remove 
the K'nim operators by commuting with finim operators, the usual orthonormalization 
condition results: 



6 Equal-time anticommutation relations 

The equal-time anticommutation relations give relations among the Haag ampli- 
tudes. These relations follow from the vanishing of the coefficients of each (linearly 
independent) normal-ordered product of in fields. Most of the relations involve Haag 
amplitudes for terms with higher-degree normal-ordered products than we have con- 
sidered, however, for the equal-time anticommutator [e,p]+ = 0, there is a relation 
that involves the wave functions: 



Thus the Haag amplitude for the off-shell electron is related simply to that for the 
off-shell proton, and the apparent asymmetry in the treatment of the constituents 
of the bound state because one particle is on-shell and one particle is off-shell is 
therefore not a true asymmetry. The two amplitudes determine each other uniquely, 
and it is convenient to define 'ipnim{j) = 4'nimi''^)- 

7 Definition of asymptotic limits 

The proper asymptotic limit is a weak operator limit that constructs an asymptotic 
field of a given mass m from the neighborhood of the mass m part of the relevant 
(product of) Lagrangian fields [T2j. The asymptotic {in or out) fields for (possibly 




(27) 



C;m(r) + V^^i^(-r) = 0. 



(28) 
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composite) particles are characterized by their rest energy E, mass m, and spin J. 
Suppressing the spin in what follows, we define the asymptotic fields associated with 
the interacting field v4(x, t) as 

(out) ^ hmt>^^^ 1 i?(x - y, t - t'; Ea, mA)A{y, t')d?y, (29) 

where the limit is the weak limit of the smeared operators, and 2^(x, t) is defined in 
Eq.(|9]). The asymptotic fields for other interacting fields are defined in an analogous 
way. (We define the asymptotic fields for composite particles in Sec. 8.) The 
asymptotic limit in momentum space is often useful in calculations: 

^in (out) ^^ ^ iimt,^^^5iE - — ) / rfE'e*(^-^')*'i(k, E') . (30) 

2m J 

Either form of the definition of the asymptotic limits makes clear that asymptotic 
fields obey the free equations of motion: 

■Q^j^tn {out)^^^^^ = {E- —V^)A"^ (''"*Hx,t), (31) 

2m 

as well as the free field anticommutation or commutation relations: 

[A*" (°"*)(x,t),A*" (°"*)(y,t')]± = ^^(x-y,t-t')- (32) 

Note that: 

r'(x,0) = 5(x) (33) 

for all E, m. 



8 Construction of asymptotic fields for bound states 

In this section, we show how to construct the asymptotic fields for the bound states 
of the hydrogen atom from a convolution of the wave function with products of 
electron and proton fields at separated points. Since the hydrogen atom is made 
of an electron and a proton and these fields anticommute, as an intuitive guess we 
start with one-half the commutator of the electron and proton fields, 

1, 
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nlm 



m + M 
mx + My' , 



3^' 



e^"(x,0)+ E / rnlrn{^-y>'''\y',t)K^lU '^'' J (34) 
n'l'm' ^ 

After normal-ordering, the terms linear in K'^'Vm' 

\Wn'Vm'{'^ - y) - V^n'Z'm'(y " x) ] /l^W ( ^ "^l'^ ) = 

n'l'm' ^ 

2^ X - y — Ti- , 35 

n'l'm' 

where we used Eq.( l28ll . We can use the orthonormality of the hydrogen wave func- 
tions to eliminate the wave functions and isolate the in field for the states of the 
hydrogen atom: The term linear in the hydrogen in field becomes: 

, * / N , / ^-, ir, / ""^x + My , ,o , , , i„ / ■'^x + My , , , 

€z™(x - y)^n'zw(x - yW:.,A — TXr)^ - = — rxr)- (36) 

m + M m + M 

Now we define an interpolating field for the hydrogen states, 

h^iU^,t) = i| d^rij^mi^MR- -^r,t),e{R+ (37) 

From these calculations, we see that the only term in this interpolating field that is 
linear in /i^"^ is just /ij^^, and further this is the only term with a singularity at the 
energy and mass of the bound state. Therefore this term can be isolated by taking 
the weak limit: 

/ili"*^(R, t) = Umr^^^ f d^R V{K -R',t- t'; E^im, m + M)hnim(R', t). (38) 

Jt' = T 

This discussion of bound states, which involves the bound-state amplitude, is an 
alternative to the constructions of bound-state operators from products constituent 
fields at a point that were formulated by Nishijima [8] and by Zimmermann [9]. 



9 Summary 



We derived the Schrodinger equation for bound states of the hydrogen atom using 
the operator equation of motion for the electron (or proton) field together with the 
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Haag expansion of the fields in terms of normal-ordered products of in (or out) 
fields. This was done without any reference to the Bethe-Salpeter equation. The 
usual normalization conditions and probability interpretation of bound state ampli- 
tudes are valid both in the present nonrelativistic example as well as for relativistic 
theories. 

The N quantum approach has been applied to many other problems, and a 
survey of results appears in Ref.( [15] )■ 



10 Appendix: the asymptotic limits 

Here we discuss the asymptotic limits of the N-quantum approach, and indicate how 
the weak limits eliminate terms that do not have delta-function or principal-value 
singularities on the energy shell of the asymptotic field. Considering an operator 
0{x,t), we define the asymptotic limits C*"*^""*) (x, t) by 

C^"(°"*)(x,t) =limt>^^^J d^x'V{^ - ^' ,t - t'; E,m)0{x' ,t'), (39) 

Using a four- dimensional Fourier transform, 

0(x,t) = J dud^kd{uj,k)e-''''+'^ '', (40) 

we find the momentum-space form of the weak asymptotic limit: 

^in{out) . ^ _ iirrn,^^^5(u -E-—) [ rfa;'e'("-"')*'0(a;', k) (41) 

2m J 

As a matrix element between two states, this gives a relation between distributions. 
If the matrix element of (9(a;',k) on the right hand side is in L^, the Riemann- 
Lebesgue lemma states that the term will vanish in this limit except when there is 
a singularity in the matrix element of 0{uj', k) on the energy shell to' = E + 
The limit is finite if the singularity is either a delta function in u' — E — ^ or has 
a principal value 1/{uj' — E — The limit diverges for stronger singularities on 
the energy shell. Therefore the requirement that the weak asymptotic limit exist 
eliminates terms that do not have delta-function or principal-value singularities on 
the energy shell. Ref.([T2]) gives further discussion of weak asymptotic limits. 
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